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We study the effect of the ferromagnetic (FM) contacts on the spin accumulation in the lateral 
spin valve system for the collinear magnetization configurations. When an additional FM electrode 
is introduced in the all-metallic lateral spin-valve system, we find that the transresistance can be 
fractionally suppressed or very weakly infiuenced depending on the position of the additional FM 
electrode, and relative magnitudes of contact resistance and the bulk resistance defined over the 
spin diffusion length. Nonlocal spin signals such as nonlocal voltage drop and leakage spin currents 
are independent of the magnetization orientation of the additional FM electrode. Even when the 
additional contact is nonmagnetic, nonlocal spin signals can be changed by the spin current leaking 
into the nonmagnetic electrode. 

PACS numbers: 72.25.-b, 73.40.Gk 



I. INTRODUCTION 

Electrons are characterized by their quantized spin 
and charge. In conventional electronic devices, only the 
charge degree of freedom has been employed for the con- 
trol of the electron transport. A new field of spintronicai 
was born from both experimental and theoretical efforts 
to harness the electron's spin degree of freedom in or- 
der to control the electric current in the devices. One of 
typical spintronic devices is the spin valve which is a hy- 
brid structure of ferromagnetic (FM) metal/nonmagnetic 
(NM) material/FM metal. The current passing through 
the spin valve depends on the magnetization configura- 
tion of two FM metals. In the collinear case, usually 
more current flows through the spin valve in the parallel 
configuration than in the antiparallel configuration. Dif- 
ference in resistance between the two is called magnetore- 
sistance. In the noncoUinear case or when two magnetiza- 
tion orientations are neither parallel nor antiparallel, the 
spin polarized current from one FM electrode exerts the 
spin torque^''^'"' on the other FM electrode, and induces 
the magnetization dynamics. Examples of a spin valve 
are giant magnetoresistance (GMR) devices^^i^ magnetic 
tunnel junctions^^ nanopillarSf^ etc. 

In contrast with vertical spin valves, lateral spin valves 
are characterized by their multi terminal functionalities 
and so are more favorable for integration into semicon- 
ductor electronics. Due to increased spacing between ter- 
minals, efficient spin injection and detection have been 
a very hot issue. The spin injection and detection ex- 
periments in the two-terminal geometry are obscured by 
other effects like anisotropic magnetoresistance. Hall ef- 
fect, etc. This defect was overcome by adopting the non- 
local spin valve geometry^ similar to the schematic device 
structure in FigU The original spin valve devices con- 
tain two FM electrodes (vs. three FM electrodes in Fig. 
[T]) contacting the nonmagnetic base electrode. In this 
lateral spin valve system, the spin transport was clearly 



observed with Al wires^" by spatially separating the spin 
current path from the charge current path and thereby 
removing other undesirable effects. The spin polarized 
current flows from the left of N (base electrode) into Fl. 
That is, spin polarized electrons are injected from Fl into 
base electrode N and is drained to the left of N. Due to 
asymmetry of two spin states in FM, the number of in- 
jected spin-up and spin-down electrons is different. In 
addition to charge current in the left of N (a; < Li), dif- 
fusion of injected spins generates spin current flowing to 
left and right of N symmetrically. Pure spin current to 
the right of N was detected^i^Siiii with another FM elec- 
trode by measuring the spin-dependent nonlocal voltage 
drop. The nonlocal spin injection and detection tech- 
nique was also used to observe^ii^ the (inverse) spin Hall 
effect in diffusive nonmagnetic metallic strips. In these 
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FIG. 1: (color online). Schematic display of the lateral spin 
valve system with three ferromagnetic electrodes, (a) Ferro- 
magnetic electrodes are labeled as Fl, F2, and F3 from left to 
right. The base electrode is denoted as N. Spins are injected 
from Fl to N by the spin-polarized current I flowing from 
the left end of N into Fl. The electrode Fl is called the spin- 
injecting probe, while Fi with 17^! will be called the nonlocal 
(voltage) probe in this paper, (b) The one-dimensional model 
geometry of the spin valve system in (a). 
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experiments, the spatial separation of charge and spin 
currents as weh as the efficient spin injection are essen- 
tial to observing the charge Hall voltage induced by the 
spin current. 

Recently experimental groups^ii^ii^iii studied the 
spin transport in the lateral spin valves with the three 
FM electrodes as shown in Fig.[T] Fl is the spin- injection 
electrode (Fgi) as usual, while F2 and F3 are the nonlocal 
voltage probes located outside the charge current path. 
While the spin polarized electrons are injected from Fl 
into N and are drained into the left of N, the transresis- 
tance is measured between Fl and F2 or between Fl and 
F3. In the former case, F2 is the spin-detecting electrode 
Fgd and F3 is the additional electrode Fa. In the latter 
case, the roles of F2 and F3 are switched. The important 
issue is how sensitively the transresistance is affected by 

Fa. 

One experimental group found that Fa can drain the 
spin current and thereby significantly suppress the tran- 
sresistance. They concluded that such additional FM 
electrode is relevant to the spin injection and accumu- 
lation in the multi terminal lateral spin valves. On the 
other hand, another groupie found that the transresis- 
tance is weakly affected by Fa even when the contact 
resistance between Fa and N is Ohmic. Moreover the 
transresistance was observecti^ to be independent of the 
magnetization orientation of Fa (parallel or antiparallel 
to that of Fsi). They concluded that such additional FM 
electrodes are irrelevant to the spin injection and accu- 
mulation in the multi terminal lateral spin valves. The 
existing experimental results seem to be contradictory to 
each other. 

Motivated by this experimental situation, we study 
theoretically the spin transport in the lateral spin valve 
with three ferromagnetic electrodes as schematically 
shown in Fig. [1] We adopt the semiclassical spin drift- 
diffusion (SDD) equation^iii^ for the one-dimensional de- 
vice structure and study the mutual effect of FM elec- 
trodes on their nonlocal spin signals such as nonlocal 
voltage drop and leakage spin current. We find that 
the efficiency of the spin current leakage into FM elec- 
trodes depends on the relative magnitude of junction re- 
sistance and the bulk resistance (defined over spin diffu- 
sion length) in FM and N electrodes. The voltage drop 
in Fsd is proportional to its leakage spin current with the 
proportionality constant given by the effective spin re- 
sistance which depends on the magnetization orientation 
of Fsd- Nonlocal spin signals are sensitive to the posi- 
tion of Fa relative to the positions of Fgi and Fgd- When 
Fa is located in between Fgi and Fgd, the transresistance 
can be either significantly or weakly affected by Fa de- 
pending on the relative magnitude of junction resistance 
and spin resistance. The effect of Fa is weak when Fa 
is located outside the region between Fgi and Fgd- Even 
though the magnitude of the spin current and nonlocal 
spin signals may be modified by Fa, the flow direction of 
the spin current in the whole device is set by the mag- 
netization orientation of Fsi and so the nonlocal spin sig- 



nals are independent of the magnetization orientation of 
Fa. This surprising result is already observedi^ in experi- 
ments and is the direct consequence of no charge current 
in Fa. Based on decoupling of charge and spin modes 
in the SDD equation and the Kirchhoff rules at the junc- 
tion, we also show that the relationship between nonlocal 
spin signals and magnetization holds true even in realis- 
tic three-dimensional samples. These interesting proper- 
ties in fact originate from zero charge current in nonlocal 
voltage electrodes. Irrelevance of magnetization orienta- 
tion of additional FM electrode to nonlocal spin signals 
implies that even additional nonmagnetic electrode can 
modify nonlocal spin signals in the spin-detecting probe. 
Our theoretical study may be useful for clarifying the 
conflicting roles^^ii^ii^iii of an additional FM electrode 
in the lateral spin valve devices. In addition, our study 
is relevant to device applications, because the multi ter- 
minal functionality is essential for device applications of 
lateral spin valves. 

The rest of this paper is organized as follows. In Sec. 
II, the spin drift-diffusion equation is briefly introduced 
and the detailed algebras for the lateral spin valve with 
three FM electrodes are included. The results of our work 
for spin valves relevant to experiments are presented in 
Sec. III. In Sec. IV, our work is summarized and its rele- 
vance to experiments is discussed. Some algebraic details 
and interesting results are included in three Appendices. 



II. FORMALISM 

From now on we are going to confine our discussion 
to the collinear magnetizations of three FM electrodes in 
spin valves and so we consider the spin polarized trans- 
port in a steady state. The noncoUinear magnetizations 
go beyond the scope of our paper since the current flow 
in the noncollinear magnetizations generates the spin 
transfer torque and induces the magnetization dynam- 
ics. In the collinear and diffusive transport, the spin 
drift-diffusion equations^iii^ have been very useful for un- 
derstanding phenomenologically the spin polarized trans- 
port in the spin valve systems. Later the SDD equations 
were derived^^ from the semiclassical Boltzmann equa- 
tion under the assumption that the spin diffusion length 
(SDL) is larger than the mean free path (MFP). Using 
the numerical solution of the spin dependent Boltzmann 
equation, the validity of the SDD equations was further 
extended'^'^ to the case when the SDL is comparable to 
the MFP. The SDD equations have been widely used for 
analyzing the spin injection experiments in various de- 
vice geometries. The SDD formalism was also applied 
to the study of spin transfer torqu o^^i^°'^^ in the case of 
noncollinear magnetizations. 

The SDD equations in the collinear magnetizations are 
written down for the spin-dependent electrochemical po- 
tential fia and electric current density Jq. Here a — ± 
represents the spin-up (-I-) and spin-down (— ) states, re- 
spectively. The presence of the spin flip scattering in 
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bulk mixes two spin states and the SDD equations can 
be written down in a matrix form. 
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Here Da is the diffusion constant for spin direction a — 

± and T-i is the average spin-flip time for an electron 

from the spin direction + to — . CTq, is the conductivity 
for electrons with spin a and e is the absolute value of 
electron charge. 

The matrix differential equation for the electrochemi- 
cal potential can be solved^'* by analyzing the eigenvalues 
and eigenvectors of the matrix in the SDD equation. One 
eigenvalue is and the corresponding eigenvector is ( J ). 
The other eigenvalue defines the spin diffusion length A 
and is given by the expression 
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and its eigenvector is 
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To find this form of the 



eigenvector, the Einstein relation as well as the detailed 
balance relation are already invoked. Obviously the first 
eigenvector (charge mode) does not discriminate between 
two spin states, while the second one (spin mode) does. 

In this section we analyze the spin polarized transport 
in the spin valve system based on the one-dimensional 
SDD equations. The device structure is displayed in 
Fig. [IKa) , where the base electrode is contacted to three 
ferromagnetic (FM) electrodes. Our primary goal is to 
understand the mutual influence of the ferromagnetic 
electrodes on the nonlocal spin signals such as the volt- 
age drops and the leakage spin currents. FM leads are 
labeled as Fi with i = 1,2,3 from left to right. The 
one-dimensional geometry, corresponding to the device 
structure, is displayed in Fig. W[h), where the junctions 
between the base electrode and the FM leads are labeled 
as X ^ Li {i = 1, 2, 3). 

In experiments, the base electrode is nonmagnetic, but 
we are going to consider the case of magnetic base elec- 
trode with its nonzero bulk spin polarization f3. Nonmag- 
netic case is recovered by a simple replacement /3 = 0. 
The spin polarization (SP) in each FM lead is denoted 
as (3i which is defined by the spin asymmetry in the spin 
dependent conductivity ai±. 



f3^ 



(4) 



With the total conductivity Ui = ai+ + crj_, the spin-up 
and spin-down conductivities can be written as 



(5) 



For the base electrode, the spin polarization (/?) in con- 
ductivity and the spin-dependent conductivities {a± ) are 
defined in a similar manner. 



When the spin polarized electrons are injected from 
Fl into the base electrode N and is drained to left, the 
electrochemical potential in the FM leads can be written 
as 
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Here Si^i is the Kronecker delta function. The spin de- 
pendent current is determined by the equation 



e dy 



(7) 



The charge current is given by the expression I^i — Ii+ + 
li- — ISi^i and fiows only in Fl but not in Fi with i ^ 
1. The i-th FM lead is contacted to the base electrode 
at X = Li. Ai, Xi, (Ti, and Pi are the cross sectional 
area, the spin diffusion length, conductivity, and bulk 
spin polarization in conductivity of the i-th FM lead, 
respectively. TZi, defined by the relation 



(TiAi 



(8) 



is the resistance of the FM electrode over the spin diffu- 
sion length. Due to an exponential decay of spin current, 
this definition of resistance makes sense physically when 
discussing the spin current. Vi is the voltage drop at 
each ferromagnetic electrode far away from the junction 
with the base electrode and is induced by the nonequilib- 
rium spin injection and diffusion. Note that the common 
Fermi energy is dropped in writing the electrochemical 
potentials in this paper, because the overall constant en- 
ergy shift does not change physics. The spin current in 



Fi I! ^ I, 



i+ 



li- is given by the expression 
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(9) 



The first term is the spin-polarized driving current, while 
the second comes from the spin accumulation and dif- 
fusion. Though no charge current flows in the region 
X > Li, the spin current is induced in the base elec- 
trode due to the spin injection, accumulation and diffu- 
sion. The spin current decays exponentially over the spin 
diffusion length and in turn leaks into the other FM elec- 
trodes, measures the magnitude of this leakage spin 
current at the interface between the base electrode and 
Fi. The leakage spin current also decays exponentially 
over the SDL in the FM electrodes. The set of six un- 
known parameters {Vi,Ii} are to be determined by the 
Kirchhoff rules at the junctions. Note that Vi and Li are 
null for 1^1, when the injected current is not spin po- 
larized. Vi can be nonzero for the tunneling barrier even 
when the injected current is not spin polarized. Hence 
we may call {Vi, li} for i ^ 1 as the nonlocal spin signals. 
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In the common base electrode, we have the electro- 
chemical potential for spin-up(+) and spin-down(— ) elec- 
trons 



e Uj.- 



J,7^ e 



Here 9{x) is the step function and the spin dependent 
current is computed from 



/+ = 



e dx 



(11) 



A, A, a, and P are the cross sectional area, the spin dif- 
fusion length, conductivity, and bulk spin polarization 
in conductivity of the base electrode, respectively. TZ, 
defined by the relation 



n 



(12) 



is the resistance over the spin diffusion length in the base 
electrode. The charge current Ic — 1+ + 1- = IO{Li — x) 
flows only at the section x < Li and the spin current 
Is = 1+ — I- is given by the equation 

If, = l3ie{Li~x)+Y,J,sgn{x-U) e-\^-^^\'\l?,) 



The additional three unknown parameters {Ji\ are in- 
troduced for the base electrode. Ji represents the re- 
distribution of spin current in the base electrode due to 
its leakage into the voltage probes (FM electrodes with 
i = 2,3). 

The electrochemical potentials are constructed such 
that the charge current is conserved at any junction in 
the device (charge continuity equation). No net charge 
current flows to the right of a; = Li and the currents 
for spin-up and spin-down electrons are equal in their 
magnitude and opposite in their sign or flowing direc- 
tion. This symmetry in spin current is strictly obeyed 
in the nonlocal spin valve geometry even in the presence 
of the multiple ferromagnetic electrodes to the right of 
the spin- injecting FM electrode Fl. Obviously the spin 
currents li and Ji are modified by the presence of other 
FM electrodes, which is our main research interest in this 
work. 

There are nine unknown parameters {Vi^Ii, Ji\ {i — 
1, 2, 3), which should be determined by the boundary con- 
ditions or Kirchhoff rules at the junctions. As already 
mentioned in the above, the electrochemical potentials 
are constructed such that the charge current is conserved. 
In our work, the spin flip scattering is neglected at the 
interface, but is taken into account in bulks. In this case, 
the spin current is conserved at each junction and the 
following relations are obtained 
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We consider the case of dirty interface between the fer- 
romagnetic electrodes and the base electrode. Due to 
a finite value of resistance at the interface, the electro- 
chemical potential across the junction is not continuous 
and should be determined by the Ohm's lawi ^^i^^i^^ 



1 ^ 

e 



/^a(y-0+)7^t: 



(15) 



Here A/iia is the difference of the electrochemical po- 
tentials at X = Li between the base electrode and the 
Fz electrode. TZti± is the spin dependent junction resis- 
tance between the base electrode and Fi, and is defined 
in terms of the spin polarization 7, of junction resistance. 
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TZti is the total junction resistance or TZu — 
TZti+Tl-ti- / (TZti+ + TZti-)- The clean or transparent con- 
tact can be recovered by a simple replacement TZti = 0. 

For the sign of /3's (spin polarization), we are going 
to adopt the following convention. When the spin-up(- 
down) electrons belong to the majority (minority) channel 
at the Fermi level, the sign of /3's is positive. On the other 
hand, the sign of /3's is negative when the spin-up(-down) 
electrons belong to the minority (majority) channel. Ac- 
cording to our convention, the sign of /3's is reversed un- 
der the magnetization reversal. The same convention ap- 
plies to the sign of 7's which are introduced to define the 
spin polarization in the resistance of the interface. 

After some algebra as detailed in Appendix El we find 
the expressions for li , Ji and Vi , which contain all the in- 
formation about the spin polarized transport in the one- 
dimensional spin valve. 



h 
I 

I 

I 



-(/3i - + Ga[(/3i - (3)Ri + (71 - !i)Rm) 
-^G.a[(/3i-/3)i?i-H(7i-/3)i?ti], (18) 

-[(/3 - (iifRi + iP^ - 2/371 + l)i?ti]5,,i 

+ [(/3,-/3)i?,; + (7,-/3)i?t,]G,i 

x[(/3i-/3)i?i + (7i-/3)i?ti]. (19) 

Gij is the element of the matrix G defined in Eq. (|A12[) of 
Appendix lAl and has the dimension of conductance. The 
other set of material parameters, Ri and Ru, are intro- 
duced in Appendix |^ and their definitions are repeated 
here for readers. 



Ri 



1-/3/ 



R = 
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These new material parameters need our special atten- 
tion. They have the dimension of resistance and deserve 
their own terminology. They are already called the spin 
resistance in the literature. First of all, the spin resis- 
tance is introduced to simplify the algebra as shown in 
Appendix [X] As the above equations show, this spin re- 
sistance determines the nonlocal spin signals such as the 
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voltage drops and the leakage spin currents in the voltage 
probes. More physical insights on the spin resistance are 
elaborated on in Appendix [b1 

Since we are interested in the nonlocal transport mea- 
surements, we focus on the leakage spin currents and volt- 
age drops in the voltage probes (Fi with i 1). The leak- 
age spin current in the nonlocal voltage probes is given 
by the expression [Eq. (fTT)) ] 

= GaiiPl - + (71 - P)Rti]I. (21) 

This relation for the nonlocal spin current suggests 
that the conductance matrix G contains all the informa- 
tion about the mutual effect of nonlocal voltage probes. 
Since the conductance Gn does not depend on the mag- 
netization configuration of the FM electrodes, the leaking 
spin current does not depend on the magnetization orien- 
tation of voltage probes, but instead depends on the mag- 
netization orientation of the spin-injecting FM electrode 
(Fl) and the base electrode (if ferromagnetic, /? 7^ 0). 
This important observation can be understood as follows. 
The flow direction of spin current (the sign of /i), in the 
base electrode as well as in the FM electrodes, is obvi- 
ously determined by the magnetization configuration in 
the spin-injecting electrode. This means the flow direc- 
tion of spin current cannot be altered by the change of 
magnetic configurations in nonlocal voltage probes. This 
is due to the fact that the nonequilibrium spin current 
is generated by the spin- injecting electrode, but not by 
nonlocal voltage probes. Another important observation 
is that the magnitude of spin current or is not modi- 
fied under the magnetization reversal of nonlocal voltage 
probes, which in fact derives from the symmetry in the 
SDD equations. This property of li derives from decou- 
pling of spin and charge modes in SDD equations as well 
as the zero charge current in nonlocal voltage probes. 
Detailed analysis can be foimd in Appendix [Cl The re- 
lation Ji = —Ii/2 simply reflects the conservation of spin 
current at the interface between the voltage probe and 
the base electrode. 

Vi defined below is ubiquitous in the expressions of 
and V^. 

V{ = [(/3i - /3)i?i + (71 - /5)i?ti]/. (22) 

We may call Vf as the spin potential which is the source 
from the spin-injecting electrode and drives the spin cur- 
rent in the spin valve. The leakage spin current in non- 
local probes can be written as 

/, = Gun'- (23) 

The spin current in the spin valve device can be expressed 
in terms of the spin potential Vf and the conductance 
matrix G as 

= f3ie{Li - x) 

-^Y.G^^^Sn{x-U)e-^^~^^^/\ (24) 

i 

If = PJ{l~e-y/^^) + [GiiVi' +pI]e-y/^\ (25) 
Zf = GaV,'e-y/^%i^l. (26) 



The voltage drop in the nonlocal voltage probes is 
given by the expression [Eq. p^ ] 

V, = m-p)R, + ij,-l3)Ru]Ga 

x[{/3i-l3)Ri + i-fi-l3)Rti]I. (27) 

Note that the voltage drop can be written as the product 
of spin current and some sort of spin resistance as 

V, = m - I3)R, + (j^ - l3)Ru]h. (28) 

Vi is the effective measure of weighted averaging the spin- 
up and spin-down electrochemical potentials in Fi and 
so depends on the magnetization configuration of Fi. In 
addition, Vi can be understood as a shift in the Fermi 
level in order to satisfy the condition of zero charge cur- 
rent in the nonlocal voltage probes. Look at Appendix 
IB] for details. Though the leakage spin current li in Fz 
is independent of the magnetization orientations of all 
nonlocal voltage probes (parallel or antiparallel to that 
of spin- injecting probe), the voltage drop Vi depends on 
the magnetization orientation of Fi, the spin-injecting 
probe and the base electrode, but not on that of other 
voltage probes. 



III. RESULTS 

In this section our discussion is confined to the spin 
valve system with three ferromagnetic electrodes and 
nonmagnetic {(3 = 0) base electrode. The current / is 
injected from the left of nonmagnetic base electrode and 
is drained to Fl (See Fig. [T|). Though the charge current 
is null to the right of the contact between Fl and the base 
electrode, the finite spin current is induced everywhere 
by the spin injection, accumulation and diffusion. With 
two nonlocal FM electrodes (labeled as F2 and F3) con- 
tacted with the nonmagnetic base electrode to the right 
of Fl, we want to study the mutual influence of two non- 
local FM electrodes on their voltage drops and leakage 
spin currents or non-local spin signals. 

With three FM electrodes, the dimension of conduc- 
tance matrix G is 3 x 3. 

/ n \Rh \Rf2\ 
G"' = iRh r2 \Rh ■ (29) 

Here r; = i?i + Ru + \R for i = 1, 2, 3 and /i = e~^^''/^, 
/2 = e-^i3/\ and h = e-^l^/^ where L„ - \L, - L.j\ 
is the distance between the contacts of the i-th and j-th 
FM electrodes with the base electrode. After inserting 
the explicit expressions of Gn into Eq. pTjl with /3 = 
or 

U ^ GaV^\ = (/3ii?i+7ii?ti)/, (30) 
we find the leakage spin currents I2 and /a in the nonlocal 
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electrodes to be given by the expressions 



IR 



i?3 + i?t3 + ^(l-e-2^=^/^)i? 



(31) 



IR 



h = -^e-^i^'/^ (/3ii?i + 7ii?*i)(i?2 + Rt2). (32) 

Here D is the determinant of the matrix G^^ and is given 
by the expression 



D = YI(r, + Ru + \r 



-i?^ ^ (.Ri + Rt 



i=l 



(-1) 



4+1 



R]ft- (33) 



Obviously the exponentially decaying factor can be ex- 



tracted out as Gil 



gn {i = 2,3), where gn is 



negative. As mentioned in the previous section, the leak- 
ing spin current {I2 and /a) does not depend on the mag- 
netization configuration of the nonlocal voltage probes 
(F2 and F3), but only on the magnetization configura- 
tion of the spin-injecting electrode (Fl). This is clearly 
explained by the fact that the spin current is generated 
by the spin-injecting electrode. 

The nonlocal voltage drops V2 and V3 are related to 
their leakage spin currents as 



V, = iP,R,+-/,Ru)h, z = 2,3. 



(34) 



This relation suggests that the nonlocal voltage drop Vi 
is proportional to the corresponding nonlocal spin cur- 
rent li and the proportionality constant is the effective 
spin resistance which depends on the magnetization ori- 
entation of Fi. This resistance is the intrinsic material 
properties of the relevant FM electrode so that the ef- 
fect of the other FM electrode is completely embedded 
into the nonlocal spin current. This means that we can 
discuss the effect of multiple FM electrodes on transre- 
sistance in terms of the leakage spin current. 

From Eqs. and ([M]), transresistance {Rs2 = 

V2/I and Rs3 — V3//) in F2 and F3 electrodes are re- 
duced to the following forms. 
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2D 
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(35) 



(36) 



We denote the transresistance in the absence of an ad- 
ditional FM electrode using the superscript as R^^ and 



R^^3 ■ The same notations with superscript will be used 
for nonlocal voltage drop and spin current. The effect 



of other nonlocal FM electrode on the transresistance 
can be quantified by computing the ratio: Rsi/R^^l'^ with 
i = 2,3. It follows from Eq. jSO]) and jMl) that 



Rs 



R 



(0) 



(0) 



AO) 



Gil 



G. 



(0) 



(37) 



The effect of an additional FM electrode on the transre- 
sistance can be measured by how much the nonlocal spin 
current is reduced or by the change of the conductance 
matrix G under the other FM electrode. 

Two important facts can be read off from Eqs. ([35]) and 
([36|) . (i) The transresistance of one FM electrode (say, 
F3) does not depend on magnetization orientation of the 
other electrode (say, F2). That is, the transresistance Rs3 
in F3 does not depend on the spin polarization, P2 and 
72, of F2. In Ref. 16, the transresistance was observed to 
be independent of the magnetization orientation (parallel 
or antiparallel to F3) of the intervening FM electrode F2, 
which is supported by our theoretical results. However 
this fact does not necessarily meani^ that the observed 
transresistance is not influenced by the additional FM 
electrode. The transresistance can be either significantly 
changed or weakly influenced by the presence of the ad- 
ditional intervening FM electrode, depending on sample 
and material parameters as we shall show below, (ii) The 
transresistance can be modified^ even when the addi- 
tional contacting electrode is nonmagnetic. The relative 
magnitude of the interface and bulk resistance (defined 
over the SDL) plays an important role in determining the 
transresistance. Irrespective of the magnetic or nonmag- 
netic nature of the intervening electrode, the transresis- 
tance will be influenced only by spin resistance and the 
interface quality. 

According to Eq. ([57]) , the effect of an additional FM 
electrode on the transresistance is equivalent to its ef- 
fect on the nonlocal leakage spin current. After the spin 
current is injected from Fl electrode, it will flow into 
both directions in N and will leak into nonlocal probes. 
From this perspective we can expect that the effect will 
be much stronger when an additional FM electrode lies in 
between two (spin-injecting and spin-detecting) FM elec- 
trodes than when it lies outside two electrodes. In the 
former case, the nonlocal spin current in spin-detecting 
probe will be reduced proportionally by the amount of 
spin current drained into an intervening electrode. In 
the latter case, the injected spin current leaks into the 
spin-detecting probe first and then into an additional FM 
electrode, so that the effect will be weaker. Mathemat- 
ically this difference between two cases comes from the 
asymmetry between Eq. ((3T|) and Eq. ([32)1 . Under the 
index exchange 2 3, /2 and I3 are inequivalent due to 
the additional term (1 - e^^^^^/^)R/2 in 

To be more quantitative, let us consider I2 when ff's 
in Eq. ([33]) are all much less than a unity. This is a good 
approximation in all-metallic lateral spin valve systems 
since the spacing between the electrodes is comparable to 
the SDL which is of the order of few hundred nanometers. 
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Under this approximation, we can readily show that 

h ^ . (38) 

Therefore, Rs2 — V2// is very weakly influenced by the 
FM electrode F3 and the transresistance is almost the 
same as that in the absence of the electrode F3. That is, 
the transresistance is not much changed by the additional 
electrode (F3) when it is contacted to the outside of Fl 
(spin current injected) and F2 (spin current detected). 
However, the effect of an additional FM electrode F3 
cannot be neglected if F2 and F3 are closer to each other 
than the SDL. So much for this case. 

We now focus on the case when an additional FM 
electrode lies in between the spin-injecting and spin- 
detecting electrodes. That is, we study the effect of F2 
on the nonlocal spin signals for F3. In the absence of the 



intervening FM electrode F2, R^^^ v. 



^25.26 



R 



(0) 
s3 



'^1^'/^ (/3ii?i -t- 7ii?ti)(/33i?3 + IzRim) 



Do 



R 
2A^' 



Ri + Rti + ^R] {R3 + Rt3 
1 



^2g-2Li3/A_ 



(40) 



leakage into F2 will be negligible and the spin current will 
mostly continue to flow in N. The leakage spin current I2 
is larger (smaller) if the effective spin resistance R2 + Rt2 
of F2 is smaller (larger), compared to the spin resistance 
R of N. Obviously the leakage into F2 reduces the spin 
current in N and in turn reduces the leakage spin current 
I3. The larger (smaller) R2 + Rt2 is, the larger (smaller) 
is I3. 

The spin diffusion length (SDL) is of the order of a few 
hundred nanometers (nm) in nonmagnetic metals and 
the SDL in FM metals is of the order of a few nm to 
a few tens of nm. The resistivity depends on the sam- 
ple quality such as the impurities, defects, etc. Though 
SDL is two orders of magnitude different between FM 
and NM, the relative magnitude of resistance {Rf- fer- 
romagnetic metal, R^: nonmagnetic metal) defined over 
the SDL can be varied from device to device. Roughly 
Rn > Rf in the spin valve devices. Usually the in- 
terface between the FM electrodes and the nonmagnetic 
base electrode is Ohmic (i?f), but not in the tunneling 
regime. In real materials, we have the following order in 
resistance: Rn > Rp > Rt- For our theoretical study, 
we will consider both cases of Ohmic and tunneling in- 
terfaces as well as other parameter regimes. 



Note also that R^^ can be obtained from Rs3 by taking 
the limit Rt2 00 or when the second intervening F2 
is effectively decoupled from the nonmagnetic base elec- 
trode. The effect of the second intervening electrode F2 
on the transresistance Rs3 can be quantified by comput- 
ing the ratio Rs3 / R^^^ , 



— [R2 + Rt2)- 



?(0) 
^s3 



(41) 



Below this general relation will be reduced to the simple 
forms case by case. 

In order to provide some physical insights, let us con- 
sider the case when ff <^\. We find the simple form of 



Rs3/R's 



R. 



s3 



i?2 + Rt 



i?2 + Rt2 



R 



(42) 



The reduction of Rs3 stems from the leakage of spin cur- 
rent into the intervening electrode F2. The efficiency of 
spin leakage into F2 is quantified by the relative magni- 
tude of the serial resistance i?2 + Rt2 in F2 and the resis- 
tance R in base electrode over the spin diffusion length. 
We can understand qualitatively the results of Eq. (142]) 
as follows. Spins are injected from Fl into the base elec- 
trode N and in turn diffuse into left and right of N. That 
is, the spin current flows in N and leaks into nonlocal 
probes. Just like charge transport, the spin current at 
the junction with F2 will continue to flow in N and also 
leak into F2. If the effective spin resistance R2 + Rt2 of 
F2 is much larger than the spin resistance R of N, the 



A. Clean F/N interface 

We consider the clean interface between the base elec- 
trode and the FM electrodes: Ri,R ^ Ru- To get the 
simple expression of Rs3, we take the limit Ru — 0. Sup- 
pose that the FM electrodes are the same material with 
roughly the same Ri ~ Rp for i — 1,2,3. If the expo- 
nentially decaying factors {ff) are negligible, we find the 
simple form of the transresistance ratio 



R 



s3 



2Rf 



R 



(0) 



s3 



2Rf + R 



R 



(0) 
s3 



-P1P3 



2i?i?p 



(2i?F + i?)2 



-ii3/A 



(43) 
(44) 



When Rf R, the transresistance will be strongly sup- 
pressed by the additional intervening FM electrode. On 
the other hand, the transresistance will be fractionally re- 
duced when i?^ is comparable to R. In the other extreme 
case oi Rf ^ R, the transresistance won't be affected by 
the intervening FM electrode. 

As noted in the above, the transresistance can be af- 
fected by the nonmagnetic electrode, /?2 = and R2 = 
Rn- Let us study this case in detail. In the clean limit 
of interface, 



R 



s3 



R 
R 



(0) 



s3 

(0) 
s3 



2Rn 
2Rf + R 

-P1P3 



2RRp 



-L13/A 



(2i?F 

Since Rn = 92^2/^2 with (32 = 0, we obtain the simi- 
lar result as in the previous paragraph depending on the 



i?)2 



(45) 



(46) 
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relative magnitude of Rn, Rp- 

If the contacts between the base electrode and Fl, F3 
are clean, but the contact with the intervening electrode 
F2 is in the tunneling regime, the effect of an additional 
electrode on the transresistance is negligible. 



B. Tunneling F/N interface 

When the junction resistance is dominant compared 
to the resistance over the spin-diffusion length in the FM 
lead and the base electrode, or when Ru 3> Rj,R, the 
expressions of the voltage drop, Eqs. |35)) and ((36)) . are 
simplified as 



^2 „ 


R 


/ 


- -^71726 




R 


/ 


- -^71736 



-Ll2/\ 



(47) 
(48) 



The voltage drop at each junction is not influenced by 
the presence of the other FM leads, when the junctions 
lie in the tunneling regime. In general, the expression 
of V3 is not affected by the presence of the second FM 
lead (additional FM lead) as far as the contact is in the 

tunneling regime. When Rt2 > Ri, Rs3 = R^s so that 
the second FM lead is effectively disconnected from the 
base electrode. 

When the accumulated spin is diffused efRciently into 
the second intervening FM lead, its effect may not be 
negligible. We still assume that the contacts with Fl and 
F3 lie in the tunneling regime. Let us see the extreme 
case of a transparent contact of F2 electrode to the base 
electrode. In this case, we may set Rt2 — and the 
desired voltage drop is given by the expression 



R 



s3 



2i?2 



R 



(0) 
s3 

?(0) 



2i?2 

R 



Rs3 



-71736 



-L13/A 



(49) 
(50) 



That is, the transresistance can be changed by the second 
intervening electrode F2 if F2 is in clean contact with the 
base electrode or if spin leakage into F2 is efficient. 



IV. DISCUSSION AND SUMMARY 

Using the one-dimensional spin drift-diffusion equa- 
tions, we studied theoretically the mutual effect of ferro- 
magnetic electrodes on non-local spin signals (the leakage 
spin currents and the voltage drops) in the lateral spin 
valve with three ferromagnetic electrodes. We found the 
generic expression of the leakage spin current, Eq. (pij) . 
and also a very simple relation, Eq. ((28)) . between the 
nonlocal voltage drop and the leakage spin current. 

Eq. ([57]) tells us that the effect of an additional elec- 
trode on the transresistance can be discussed in terms of 



the leakage spin current and in turn in terms of the con- 
ductance matrix. The measured non-local spin signals 
depend on the position of an additional FM electrode rel- 
ative to the spin-injecting and spin-detecting electrodes. 
When the additional electrode lies outside of the two FM 
electrodes, non-local spin signals are found to be weakly 
influenced due to the exponentially decaying spin coher- 
ence over the SDL. On the other hand, when it is lo- 
cated in between the two FM electrodes, the non-local 
spin signals can be strongly modified provided the junc- 
tion resistance is lower than or comparable to the spin 
resistance defined over the spin diffusion length in the 
FM electrodes and the nonmagnetic base electrode. If 
the junction resistance is high, the non-local spin signals 
are weakly modified even when the additional FM elec- 
trode is located in between the two FM electrodes The 
most general expression for the transresistance ratio is 
given by Eq. (j4ip . In general, the non-local spin signal 
is not much modified when the additional electrode is in 
tunneling contact with the base electrode, but is frac- 
tionally reduced when the contact is Ohmic. We also 
found that the non-local spin signals are independent of 
the magnetization orientation of the additional FM elec- 
trode, which agrees with the experimental observationi^^ 
This result suggests that even the intervening nonmag- 
netic electrode can change non-local spin signals, which 
was already observed^^ experimentally. 

Since our theoretical study is based on the one- 
dimensional device structure, some care is needed when 
we try to apply our theoretical results to interpreta- 
tion of experimental data. Strictly speaking, the ex- 
perimental spin valve structure is not one-dimensional 
in terms of the current distribution. Harmle et al. nu- 
merically showed^! that the nonlocal voltage drop de- 
pends strongly on the spatial distribution of the spin- 
polarized current. The one-dimensional approximation is 
valid when the current is uniformly distributed through 
the contact. When the contact is clean between the FM 
electrode and the base nonmagnetic electrode, the cur- 
rent flow may well not be uniform through the interface— 
and may be short circuited. In this case, the nonlocal spin 
signals may deviate from its theoretical estimate based 
on one-dimensional SDD equations. Keeping these re- 
strictions in mind, let us apply our theoretical results to 
two experimental worksJ^^^^ 

For numerical estimation {Rs/Rso = Rss/R'^^^ in this 
section), we take examples of Co/Cu/Co and Py/Cu/Py 
lateral spin valves and use the following sample size and 
material parameters. The thickness and width of the 
nonmagnetic base electrode are taken as 80 nm and 300 
nm, respectively. The width of all the ferromagnetic lay- 
ers is assumed to be the same as 100 nm. The sepa- 
ration between nearest ferromagnetic layers is taken as 
200 nm, which gives 300 nm of center-to-center distance. 
We use material parameters measured at low tempera- 
tures. The parameters for Cu are l/ccu = 6x 10^^ Jim — 
and Acu = 1 /xmiiS For Co, we use Pco = 0.46,— 
7CO/CU = 0.77,5^ l/aco(l - f3^co) = 7.5 x 10-« Qm^ 



9 



Aco = 59 nm^ and i?co/Cu^/(l - 7Co/Cu)^ = 0.52 x 
10-15 nm^^ For Py, we take /3py = 0.73^ 7Py/Cu = 
0.70,^ l/crpy(l - /3py) = 15.9 X 10-8 flm^ Apy = 5.5 
nm,^and i?Py/cu^/(l -7Py/Cu)' = 0.54 xlO-^^ Om^^ 

For the Co/Cu/Co spin valve, R/2 = 125 mil, R2 = 
150 mil, and Rt2 — 17 mil are obtained. The estimated 
spin signal is reduced to the value Rg/Rso = 0.68 by 
the intervening F2 electrode. V2/I is also reduced by 
the factor 0.87 due to F3 electrode. For the Py/Cu/Py 
spin valve, we have R/2 = 125 mil, R2 — 29 mil, and 
Rt2 = 18 mil. The reduced spin signal V3// by the F2 
electrode is Rs/Rso — 0.46. V2// is reduced by the factor 
0.87 due to F3 electrode. 

Since the SDL of Co is rather long, R2 is comparable 
to R/2 in the Co/Cu/Co case and Rs/RsO is large. Since, 
in the Py/Cu, R/2 is bigger than R2 and Rt2, Rg/Rso 
is small. The rather significant reduction of estimated 
V2/I in both cases stems from our choice of the long 
SDL of Cu at low temperatures. The long SDL means 
that the chemical potential splitting between opposite 
spin directions, though exponentially decaying, remains 
significant up to the position of the F3 electrode. The 
significant leakage of spin currents into F3 results in re- 
duction of the spin signal. At room temperature, the 
SDL of the base electrode (nonmagnetic metal) is a few 
hundreds nanometers such that the reduction of V2 // by 
the F3 electrode is only a few percents. For the experi- 
mental conditions in Refs. [3 and [161 R/2 is comparable 
to i?2 + Rt2 and we can estimate theoretical value of 
Rs/RsO, Eq. dUD: Rs/RsO - 0.5 although the observed 
Rs/Rso is smaller for Ref. [3 and is close to a unity for 
Ref. E 

As pointed out in Ref. [l^ the contact between the 
Permalloy electrode and the base Ag wire is very clean 
and the point injection and detection of current is 
suggested. In this case, the current distribution in 
the devices may well be nonuniform so that our one- 
dimensional theory cannot be straightforwardly applied. 
We believe that the nonuniform current distribution is 
the main reason why some of our theoretical estimates 
are in poor agreement with the results of Ref. We 
may discuss the relevance of the device dimensionality 
based on the effective spin resistance. According to our 
theoretical analysis, nonlocal spin signals are determined 
by the relative magnitude of junction resistance and spin 
resistance in FM and NM electrodes. This relevant resis- 
tance is defined under the assumption that the current 
distribution is uniform in the device. When the current 
distribution is not uniform as in real devices, we may still 
be able to define the spin resistance using the effective 
cross sectional area which is smaller than the geometrical 
cross section of the sample. Nonuniform current distri- 
bution tends to increase junction resistance as well as 
spin resistance, and will modify the magnitude of nonlo- 
cal spin signals. This may be one of the reasons for the 
discrepancy between two experimental results J^ii^ 

Nonlocal spin signals (the leakage spin current and the 
voltage drop) in one nonlocal FM electrode are shown 



not to depend on the magnetization orientation (parallel 
or antiparallel) of the other nonlocal FM electrode. We 
believe this symmetry of nonlocal spin signals are robust 
against the sample dimensionality, though their magni- 
tude is sensitive to samples. The spin current is generated 
by the spin-injecting FM electrode and so its flow direc- 
tion cannot be changed by the magnetization orientation 
of nonlocal FM electrodes. In addition, the magnitude 
of spin current does not depend on the magnetization 
orientation of nonlocal FM electrodes. This property de- 
rives from both decoupling of spin and charge modes in 
the SDD equations and zero charge current in nonlocal 
voltage probes. Hence our conclusion about the relation- 
ship between nonlocal spin signals and magnetization in 
nonlocal voltage probes won't depend on the sample di- 
mensionality and qualities. This point is demonstrated 
more explicitly in Appendix [Cl 

Finally, we would like to discuss the properties of tran- 
sresistance (TR) and (longitudinal) magnetoresistance 
(MR) in spin valves under magnetization reversal. Ob- 
viously, both TR and MR are modified under magne- 
tization reversal of two probing FM electrodes. Under 
magnetization reversal, TR changes its sign while MR 
changes its value. Note that MR, in general, consists 
of the two contributions: one part (background) remains 
the same but the other changes its sign under magnetiza- 
tion reversal. Let us consider the effect of an additional 
FM electrode (Fa) on TR and MR. For the vertical spin 
valves, it won't be easy to implement Fa. So we consider 
TR and MR in the lateral spin valves with Fa. Usually 
MR is obscured by other effects in the lateral spin valves 
as mentioned before. However, with increased SDL, MR 
was successfully measure d^"'^^ in the carbon nanotube 
and graphenes. To measure MR in the lateral spin valves 
of Fig. [U Fl and F3 are both current and voltage probes. 
Based on the results of Appendix [Cl we can argue that 
MR should be independent of the magnetization orien- 
tation of F2,'^^ because there is no charge current in F2 
(an additional electrode). Explicit calculation;^ using 
the SDD equations, confirms this claim. That is, both 
TR and MR are independent of the magnetization orien- 
tation of Fa. On the other hand, TR and MR depend on 
the magnetization orientation of ferromagnetic electrodes 
through which the charge current flows. 
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APPENDIX A: DETAILS IN ALGEBRA 

In this Appendix we show that the algebraic manipu- 
lation can be highly simplified by a proper definition of 
material parameters and by the vector and matrix no- 
tations. Kirchhoff rules lead to the constraints given by 
Eqs. and Eq. can be written down explic- 
itly leading to the following six relations {i = 1, 2, 3) 



,-\L.-L,\/\ _ 



= [(l±/3i)/4i±/d 



1±7. 



(Al) 



For algebraic convenience, we introduce new material pa- 
rameters 



7^^ 



R = 



n 



1-13- 



Rti — 



Tit 



(A2) 



andAij = ^ These new material parameters 

highly simplify the complicated algebra as well as deter- 
mine the spin currents in the nonlocal voltage probes. 
Then the voltage drop Vi can be written as 

j 

-(1 T l^)Ru [(1 ± f3l)I5^,l ± /,] . (A3) 
Addition and difference of two ViS lead to 

-PRY^A^jJj, (A4) 



and 

{R^ + Rti)h-RY,A^jJo = ill - f3i)RtiI S,^iA5) 

It is much more convenient to introduce the matrix no- 
tation for the algebraic manipulation. |/) = {I1I2I3Y, 
\J) = {J1J2J3Y, \V) = {ViV2V3y, and |1) = (100)*. Here 
the superscript t represents the transpose of row vectors 
so that its effect is to change them into column vectors. 
With these notations, Kirchhoff rules can be written in 
compact forms as 



\V) = 

\J) = 
= 



-(1 - /3i7i)i?ti/ |1) + [/3R + 7Rt] |/) 
-PRA\J), (A6) 

/3i)/|l), (A7) 
(R + Rt)|/) -i?A|J) 

+ {l3i-li)RtiI (A8) 



Here R and R* are diagonal matrices with diagonal ele- 
ments Ri and Ru, respectively. Similarly, /3 and 7 are di- 
agonal matrices with diagonal elements representing the 



spin polarization of each FM electrode and the junction 
resistance, respectively. A is the matrix with its elements 
given by Aij . Formally, the unknown parameters can be 
written in a more compact matrix form as 



(/3-/3i)/|l) 
+ [(/3i-/3)i?i 



1^) 



4[(A 



(71 -/3)i?ti]/G|l), (A9) 

f3)Ri + i-fi~P)Rti]IG\l), (AlO) 

[(/3 - pifRi + (/?" - 2/371 + l)Rti]m 
+ [(/3i-/3)i?i-f(7i-/3)i?ti] 
x[(/3-/3)R + (7-/3)Rt]/G|l). (All) 

Here the matrix G is defined by the expression 

-1 



G = 



R 



Rf + -RA 
2 



(A12) 



The matrix G, with the dimension of conductance, is 
independent of magnetization configurations (parallel or 
antiparallel to the spin-injecting electrode Fl) of the FM 
electrodes. The set of parameters, li, Ji and Vi, contains 
all the information about the spin polarized transport in 
nonlocal spin valves. In components, we have the spin 
currents 



i = (/3-/?i)5.,i 



I 



l3)Ri + 
~f3)Ri 



(71-mti], 
f (71 - P)Rti] 



(A13) 
(A14) 



and the voltage drops in the FM electrodes 

y = -[(/3-/3i)2i?i + (/?2_ 2/371 + l)i?ti]<5.,l 

+ m-f])R, + ij,~P)Ru] 
xG,i[(/3i-/3)i?i + (7i-/3)^ti]- (A15) 

Note that the final results are written down in a very 
compact form, using new material parameters as well as 
the conductance matrix. 



APPENDIX B: PHYSICAL MEANING OF SPIN 
RESISTANCE 

Spin resistance was defined in order to simplify the al- 
gebra. In this section we are going to infuse some physical 
meaning into spin resistance. Let us recast Eq. (|Aip for 
i =/= I (nonlocal voltage probes) into a more illuminating 
form as 



± 



l±/3 



±- 



h ( 2nt 



27^, 



1 ± 7, 1 ± /3, 



, (Bl) 



where Ui acts as the effective electric potential of the base 
electrode at the junction with Fi and is defined by 



(B2) 
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For electrons with negative charge, Vi is the electric po- 
tential for both spin directions far into Fi, and ±[7^/(1 ± 
/?) is the electric potential for spin-up and spin-down elec- 
trons, respectively, of the base electrode at the junction 
with Fi. Refer to Eqs. ([6]) and pO]) . We can deduce that 
the current Ii/2 at the interface flows into (out of) Fi for 
spin-up (spin-down) electrons. 

The left hand side (LHS) of Eq. ((BT|) represents the 
electric potential difference for both spin directions be- 
tween the base electrode and Fi at the deep inside. The 
right hand side (RHS) is the product of the current Ii/2 
and the effective spin-dependent resistance. The sign 
in front represents correctly the flowing direction of the 
spin-up and spin-down current, respectively. The first 
term in the parenthesis is the spin-dependent tunnel re- 
sistance as defined in Eq. (jl6p . The second term is none 
other than the spin resistance, which was introduced in 
the main text. With this spin resistance, Eq. (|Bip is 
the effective Ohm's law for the leakage spin-up and spin- 
down currents. 

When Ui is eliminated from Eq. (|Bip . the relation be- 
tween the nonlocal voltage drop Vi and the leakage spin 
current li or Eq. (^5]) is obtained. From Eq. (|B1[) . we 
can deduce the physical meaning of Vi. No charge cur- 
rent flows in the nonlocal FM electrodes. Vi represents 
the shift of the electrochemical potential in Fi to satisfy 
the constraint of no charge current flow. If we eliminate 
Vi from Eq. (|Bip . we find the following relation 



li {Rti + Ri 



2 ^ 



i?^ Jje-l^'-^^l/\ (B3) 



The material parameters are defined in Appendix [X] 
How do we interpret this relation? This relation can 
be considered as the Ohm's law for the leakage spin cur- 
rent. The LHS is the product of the spin current li and 
the effective resistance. From the standpoint of Fi, the 
spin current li flows from both sides of the base electrode 
[R) through the junction {Ru) and into Fi {Ri). Hence 
the effective resistance is R/2 -\- Ru + Ri as in the above 
equation. The RHS is the effective spin potential which 
combines the source term from Fl and the sink terms 
from other nonlocal FM electrodes. 



APPENDIX C: DEPENDENCE ON 
MAGNETIZATION DIRECTIONS: 
THREE-DIMENSIONAL CASE 

One of key results of our paper is the independence of 
the spin accumulation, the spin current, and the nonlo- 
cal voltage on the magnetization directions of electrodes. 
This appendix is aimed to provide an insight into the ori- 
gin of this independence in three-dimensional situations. 
We again consider the geometry in Fig. 1. Similar nota- 
tions will be used. The three-dimensional SDD equation 
is given by Eq. ([1]) and the associated charge and spin 



current densities j'^, j'^ are given by 

f = -V(a+/i+ + a_/i_), (CI) 
e 

f = iv(a+M+-^-M-)- (C2) 

e 

Similar relations hold for fiia, j^, and in the Fz elec- 
trode. The system is subject to the following boundary 
conditions. From the condition of no leakage current to 
air or insulating substrate. 



n • V/^a = fij • V/ija = 



(C3) 



should hold at the sample boundaries facing air or insu- 
lating substrate. Here n and rii denote normal vectors 
perpendicular to the boundaries. From the constraint of 
the current continuity applied to the interface between 
the base electrode and the electrode Fi, one finds that 
the following relation should hold at the interface. 



(C4) 



The Ohm's law provides another boundary condition for 
the interface. 



1 



(C5) 



When combined with these boundary conditions, the 
SDD equation completely fixes the spin-dependent elec- 
trochemical potentials. Here we remark that the current 
density VlaiaHia/e), instead of the current lia, appears 
in Eq. (|C5[) and thus TZua hi Eq. (|C5|) amounts to the 
spin-dependent junction resistance per unit area, instead 
of the junction resistance. We also remark that the tun- 
neling barrier at a junction may not be uniform in re- 
alistic experimental situations and such non-uniformity 
can be taken into account by simply regarding TZua as a 
position-dependent quantity since Eq. (jC5[) remains still 
valid even for the nonuniform barrier as long as the tun- 
nelling current remains perpendicular to the interface. 

In order to examine the dependence of the spin accu- 
mulation on the magnetization directions, we reexpress 
the involved equations in terms of the spin accumulation 
/i_ and the charge potential {a+iJ,^ + a-iJ-)/a. The 
SDD equation JT]) is again decomposed into the following 
two decoupled equations (spin mode and charge mode) (21 

V2(/i+-/i_)^l(^+-^_), (C6) 

V2(a+M+ = 0. (C7) 

The boundary conditions for the spin accumulation can 
be derived from Eq. (|C3[) and one obtains 

n • V(^+ - ^_) = fii • V(^i+ - /ii_) = 0, (C8) 

from Eq. (|C4p . one obtains 

(C9) 



n • V{ni+ - Hi-) 
a 1-13^ 



2 

l + f3 



a, 

e 



1 



1 



I + (3, 1 - (3, 
l-(3 



n- V(/i+ -/i_) 



f3^ 
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and from Eq. (|C5|) . one obtains 



from Eq. (jC4|) . one obtains 



1 



a,TCt,(l - 
4e 



n • V(/ii 



(CIO) 

+ n-j^ 



Now we are ready to discuss the magnetization direc- 
tion dependence of the spin accumulation, which is com- 
pletely fixed from its SDD equation (jC6|) and boundary 
conditions (|C8)) . ((C9ll . and (jClOp . Note that in these 
equations, all terms that depend on the magnetization 
directions are multiplied by the charge current density. 
Thus the spin accumulation should be independent of 
the magnetization directions of electrodes in which the 
charge current density vanishes. 

Next we discuss the magnetization direction depen- 
dence of the spin current. The spin current can be ob- 
tained from the spin accumulation as follows. 



J. = 



2e 
2e 



(Cll) 

V(/i.+ -A^.-) + /3j^ (C12) 



Then from the properties of the spin accumulation, it 
is evident that the spin current density should be inde- 
pendent of the magnetization directions of electrodes in 
which the charge current density vanishes. 

Finally we discuss the charge potential, which is sub- 
ject to the SDD equation (jC7p . The boundary conditions 
for the charge potential can be derived from Eqs. (|C3p . 
and one obtains 



(C13) 



n • V((T+^+ + (T-fJ.-) = n • V(o-i+^i+ 



and from Eq. (jCSp . one obtains 



-M.-), (C14) 



e CT; 



(C15) 



2e 

fT,7^t.(l-/32) 



[(i + A)(i + 7.)-(i-A)(i-7»)] 



16e 

xn • V(/ii+ - m-) 
+ ^ [(1 + + 7.) + (1 - - 7.)] A • 



Note that the SDD equation (|C7p and the boundary con- 
ditions (|C13p . (jC14p are not dependent on the magnetiza- 
tion directions of any electrodes. Thus the magnetization 
direction dependence can arise only from the boundary 
condition (|C15p . From Eq. (|C15p and from the properties 
of the spin accumulation, one then finds that the charge 
potential at the electrode Fi is independent of the mag- 
netization direction of other non-current-carrying elec- 
trodes. This in turn implies that the nonlocal voltage 
measured between the electrode Fz and the base elec- 
trode {x = -\-oo) should be independent of the magneti- 
zation directions of other non-current-carrying electrodes 
Fj (i ^ i). 
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